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Summary 

The  objective  of  this  NBS  program  is  to  identify  and,  insofar  as  practical, 
remove  technical  barriers  to  the  development  of  instrumentation  for  use  in 
detecting  and  locating  incipient  faults  in  underground  power  transmission 
systems. 

Underground  transmission  of  electrical  energy  has  economic  a"d  environmental 
advantages  compared  to  overhead  transmission.  A major  concern,  however,  is 
reliability.  Failures  are  more  costly  to  locate  and  repair. 

The  Department  of  Energy  has  a program  to  develop  instrumentation  for 
the  detection  and  location  of  incipient  faults  (the  condition  within  a cable 
n which  insulation  degradation  is  occurring  at  a significantly  higher  than 
normal  rate  which,  if  not  corrected,  will  cause  the  cable  to  fail  prematurely) . 

If  this  program  is  successful,  cables  about  to  fail  can  be  repaired  during 
scheduled  outages  minimizing  perturbation  to  the  transmission  system. 

Three  contracts  were  awarded  by  DoE  to  develop  such  instrumentation. 

Besides  monitoring  the  technical  progress  under  these  contracts,  NBS  has 
provided  technical  support  where  needed.  Technical  barriers  include  lack 
of  knowledge  of  (1)  the  physical,  chemical,  and  electromagnetic  behavior  of 
cable  insulation  immediately  preceding  breakdown;  (2)  the  manner  in  which 
characteristic  information  concerning  this  behavior  is  propagated  through  the 
power  system;  and  (3)  appropriate  methods  of  coupling  detection  systems  to 
operating  transmission  cables. 

A measurement  system  has  been  developed  at  MBS  to  enable  studies  of 
the  rf  properties  of  cables,  and  of  the  spectral  content  of  discharges  at 
incipient  fault  sites,  and  to  perform  time  and  frequency  domain  refl ectometry 
and  other  related  measurements. 


Preliminary  measurements  were  aimed  at  determining  the  sensitivity  of 
cable  sounding  techniques  (frequency  or  time  domain)  in  detecting  impedance 
changes  in  the  cable.  Changes  of  the  order  of  a few  parts  per  million  (ppm)  were 
detectable  in  some  cases.  A sharp  needle  was  inserted  into  the  cable  dielectric 
and  voltage  applied.  After  several  breakdowns  the  resulting  breakdown  channels 
were  not  detectable.  The  high  frequency  signals  necessary  to  interact  with 
small  damage  sites  were  apparently  not  reaching  the  breakdown  region.  Substantia' 
effort  was  applied  to  improving  the  coupling  between  the  TDR  (Time  Domain 
Reflectometry)  unit  and  the  cable  under  test.  While  this  improved  matters, 
the  intrinsic  attenuation  in  the  dielectric  was  shown  to  rule  out  sounding 
approaches  in  the  cable  selected  for  these  measurements.  Such  techniques 
might  still  be  useful  in  other  types  of  cables  and  for  fault  location  in  general. 

Future  measurements  will  characterize  the  rf  properties  of  underground 
transmission  cables.  Extensive  software  has  been  developed  to  permit  these 
measurements  in  the  time  domain  and  analysis  in  the  frequency  domain. 


Abstract 


Technical  barriers  exist  in  the  development  of  instrumentation  to  detect 
and  locate  incipient  faults  in  underground  transmission  cables.  Knowledge 
is  required  of  the  physical,  chemical,  and  electromagnetic  properties  of  cables 
which  precede  breakdown,  of  the  manner  in  which  characteristic  rf  signals 
propagate  in  cables,  and  of  appropriate  methods  of  coupling  detection  systems 
to  operating  transmission  cables.  A measurement  program  has  been  initiated 
that  will  provide  data  on  the  rf  properties  of  cables  and  on  the  characteristics 
of  some  forms  of  incipient  faults.  Preliminary  measurements  demonstrate  the 
limitations  of  frequency-  and  time-domain-reflectometry  techniques 
in  the  detection  of  incipient  faults.  Software  is  presented  which  permits 
the  Fourier  transform  of  step-like  waveforms. 


I.  Introduction 


The  objective  of  this  NBS  program  is  to  identify  and,  insofar  as 
practical,  remove  technical  barriers  to  the  development  of  instrumentation 
for  use  in  detecting  and  locating  incipient  faults  in  underground  power 
transmission  systems. 

In  many  situations  underground  power  transmission  cables  have  economic 
(and  environmental)  advantages  over  overhead  transmission  lines  resulting 
in  an  increase  in  underground  installations.  Underground  transmission  does 
have  a major  disadvantage  in  the  area  of  reliability.  Faults  can  be  extremely 
difficult  to  locate  accurately  and  repair  with  the  results  that  either  long 
unscheduled  outages  or  extensive  redundant  cables  are  necessary. 

To  alleviate  this  problem,  the  Department  of  Energy  has  a program  area 
entitled  "Research  and  'Development  of  Fault  Detection/Location  Techniques 
and/or  Equipment  for  Underground  Power  Transmission  Cable  Systems."  The 
output  of  this  program  is  expected  to  be  instrumentation  and  techniques  for  the 
detection  and  location  of  incipient  faults  in  different  types  of  underground 
cables;  viz.  pipe-type  (high  or  medium  pressure  oil  or  gas  impregnated), 
extruded  solid  dielectric,  compressed  gas,  and  possibly  cryogenic.  An 
"incipient  fault"  is  defined  to  be  the  condition  within  a cable  system  in 
which  local  insulation  degradation  is  occurring  at  a significantly  higher 
than  normal  rate  which,  if  not  corrected,  will  cause  the  cable  system  to 
fail  prematurely.  A reliable  incipient  fault  detector/locator  would  allow 
the  defective  area  to  be  repaired  during  a scheduled  outage,  minimizing  the 
perturbation  to  the  transmission  system. 
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The  Department  of  Energy  in  the  summer  of  1978  awarded  contracts 
to  Purdue  University,  SRI  International,  and  Westinghouse  Corporation  to 
develop  (or  determine  the  feasibility  of)  instrumentation  for  the  detection/ 
location  of  incipient  faults.  The  School  of  Electrical  Engineering,  Purdue 
University  proposed  to  develop  instrumentation  that  would  correlate,  using 
variable  time  delay,  the  noise  signal  received  at  both  ends  of  the  cable 
emanating  from  partial  discharges  at  incipient  fault  sites.  The  Remote 
Measurements  Laboratory  of  SRI  International  proposed  instrumentation 
based  on  swept-frequency  reflectometry.  A two-frequency  mixing  (resulting 
from  the  nonlinear  incipient  fault  region)  approach  was  proposed  by  the 
Electromagnetics  Sciences  Laboratory  of  SRI  International.  The  Chemical 
Sciences  Division  of  the  Westinghouse  R&D  Center  proposed  to  do  a feasibility 
study  on  the  use  of  acoustic  waveguides  for  sensing  incipient  faults.  As 
their  work  evolved  this  also  included  acousto-optic  techniques  using  fiber  optics. 

The  MBS  contribution  to  this  program,  besides  monitoring  the  technical 
progress  of  the  three  contractors  for  DoE,  was  to  provide  technical  support 
to  the  research,  development,  and  testing  of  incipient  fault  detection/location 
instrumentation.  In  general,  the  MBS  research  program  was  directed  towards 
providing  knowledge  of  the  physical,  chemical,  and  electromagnetic  behavior 
which  immediately  precede  breakdown  of  cables,  of  the  manner  in  which 
characteristic  information  concerning  these  events  is  propagated  through 
the  power  system,  and  of  appropriate  methods  of  coupling  detection  systems  to 
operating  transmission  cables. 
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II.  Technical  Progress 

II. 1 Measurement  System 

During  this  year  a measurement  system  has  been  developed  at  NBS  to 
enable  studies  of  the  rf  properties  of  cables  and  of  the  spectral  content  of 
discharges  at  incipient  fault  sites,  to  perform  time-  and  frequency-domain 
reflectometry  and  to  conduct  other  related  measurements.  Time-domain- 
reflectometry  (TDR)  measurements  were  performed  on  cables  with  induced 
defects  to  determine  the  sensitivity  of  TDR  techniques  for  incipient  fault 
detection.  Software  was  developed  that  will  permit  the  determination  of 
the  velocity  and  attenuation  of  rf  signals  in  cables  as  a function  of 
frequency  (or  equivalently  the  transfer  function)  in  the  time  domain. 

Figure  1 shows  a block  diagram  of  the  measurement  system  for  generating 
signals,  capturing  signals,  and  providing  signal  analysis  capabilities  in 
both  the  time  and  frequency  domains J The  heart  of  the  system  is  the 
Tektronix^  4051  intelligent  terminal  with  its  peripherals.  The  intelligent 
terminal  is  connected  to  both  the  RS-232  and  general  purpose  interface  bus 
(GPIB).  This  provided  direct  communication  between  other  system  components 
shown.  Also  through  the  modem  the  terminal  can  communicate  with  other 
intelligent  sources  such  as  the  NBS  Univac^  1108  computer  or  with  a 
minicomputer  at  relatively  high  data  transfer  rates.  The  terminal  outputs 
can  go  to  the  Univac^  1108,  via  the  RS-232  bus,  the  printer  via  the  RS-232 
bus,  the  plotter  via  the  GPIB  bus,  or  to  a hard  copy  unit.  Directly  attached 
to  this  system  is  the  floppy  disc  file  manager  system  capable  of  storing 
approximately  1900  kilobytes  of  information  for  direct  access.  Both  data 
files  and  programs  can  be  accessed  rapidly.  The  ROM  (read-only-memory) 
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expander  provides  eight  slots  for  specially  designed  and  programmed  "ROM 
packs"  for  editing,  matrix  functions,  binary  conversions  and  loading,  special 
signal  processing  routines,  etc.  This  adds  a quick  utility  to  the  system 
without  using  memory  overhead  or  requiring  additional  programming  time. 

Having  32  kilobytes  of  memory,  the  4051  system  possesses  substantial  "stand 
alone"  analysis  and  processing  power.  Large  jobs  are  more  efficiently 
managed  by  batch  processing  by  the  Univac^  1108  or  other  large  system.  A 
Tektronix^  4054  intelligent  terminal  is  also  available  which  is  an  order 
of  magnitude  faster  than  the  4051. 

The  Norland^  2001 A digital  oscilloscope  is  a digitizing  oscilloscope  with 
sample  times  as  short  as  2 ps/point  for  a 10-bit  word  to  1 ys/point  for  an 
8-bit  word.  It  can  digitize  simultaneously  up  to  4 input  channels  each  of 
1024  points,  2 simultaneous  channels  of  2048  points  each,  or  1 large  continuous 
array  of  4096  points.  It  has  signal  processing  intelligence  that  can  do, 
among  other  things,  FFT  (fast  Fourier  transform)  and  correlations  with  single 
keystroke  commands.  Programs  and  data  can  be  stored  and  retrieved  from  its 
floppy  disc  system.  It  currently  communicates  with  the  remaining  system 
over  the  RS-232  bus  but  can  be  adapted  to  use  the  GPIB  structure  through  an 
interface. 

The  Tektronix^  7704A  oscilloscope  coupled  with  the  Tektronix1  P7001 
Prograrroner  forms  a versatile  digital  processing  oscilloscope  (DP0)  which 
can  accept  spectrum  analyzer  plug-ins,  time-domain  reflectometer  plug-ins 
as  well  as  numerous  Tektronix^  series  7000  plug-ins.  Using  on-hand  spectrum 
analyzers,  capability  exists  for  spectral  analysis  over  a frequency  range 
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from  20  Hz  to  1.8  GHz.  Because  the  DPO  is  interfaced  onto  the  GPIB  structure, 
the  digitized  output  from  any  of  the  plug-ins  can  be  "captured"  and  processed 
by  the  intelligent  terminal.  The  digitized  waveform  consist  of  512  10-bit 
words.  A high-speed-average  option  is  included  which  permits  up  to  4096 
different  waveform  samples  to  be  averaged  before  transmission  to  the 
intelligent  terminal.  In  principle  these  4096  waveforms  could  be  transmitted 
sequentially  to  the  terminal  which  could  do  the  averaging,  but  because  the 
transfer  time  is  relatively  long,  the  latter  method  is  not  normally 
practical . 

The  Tektronix^  7912AD  is  a very  high  speed  state-of-the-art  digitizer 
capable  of  digitizing  rates  up  to  1 GHz.  Fast  transients  can  be  sampled 
every  nanosecond  for  an  array  of  512  points  and  with  a resolution  of  9-bits. 
Using  a special  direct  access  plug-in,  rise-times  as  fast  as  350  ps  can  be 
digitized.  Data  can  be  output  over  the  GPIB  to  the  intelligent  terminal. 

The  time-domain  reflectometry  equipment  consists  of  the  Tektronix^  7S12 

plug-in  and  associated  sampling,  pulse  generator,  and  trigger  recognizer  heads. 

Used  for  TDR,  the  7S12  can  measure  reflection  coefficients,  p,  from 
-3  -'3 

2 x 10  /division  to  500  x 10  "/division.  (When  p = ±1 , there  is  total 
reflection,  when  p = 0,  there  is  no  reflection).  Used  as  a general-purpose 
sampler,  repetitive  waveforms  up  to  18  GHz  can  be  displayed.  This  TDR/sampler 
can  be  used  as  a plug-in  in  the  DPO  permitting  digitizing  and  transmitting 
of  waveforms  to  the  intelligent  terminal. 

Three  spectrum  analyzers  are  available.  Two  of  them  (Tektronix^  7L5  and 
7L13)  can  be  used  in  the  DPO.  The  7L5  has  a frequency  range  of  20  Hz  to 
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5 MHz  and  the  7L13  has  a 1 kHz  to  1.8  GHz  range.  A third  spectrum  analyzer 
(Tektronix^  5L4N)  is  available  for  low  frequency  measurements  (20  Hz  to  100  kHz). 

Besides  the  above-described  instrumentation  the  measurement  system 
consists  of  various  pulse  generators,  counters,  and  sweep-frequency  oscillators 
most  of  which  can  communicate  with  the  intelligent  terminal  over  the  GPIB. 

II .2  Experimental  Results 

Time-domain-reflectometry  measurements  were  made  on  cable  samples 
containing  known  defects  in  order  to  determine  the  limitations  of  TDR  methods 
for  incipient  fault  detection.  The  TDR  unit  was  used  as  a plug-in  in  the 
digital  processing  oscilloscope  permitting  communication  with  the  intelligent 
terminal.  Software  (see  Appendix  A)  was  written  to  process  this  data.  The 
software  permits  automatic  acquisition  of  waveforms,  high-speed  averaging, 
proper  scaling,  plotting,  and  filtering.  Portions  of  the  waveforms  can 
even  be  enlarged  for  more  detailed  plots. 

A series  of  TDR  measurements  was  performed  on  a short  length  of  distribution 
cable.  The  cable  studied  was  a 6.2  m long,  2 cm  diameter,  15-kV,  extruded- 
polyethelene  distribution  cable  with  semiconducting  screen  layers.  The  cable 
had  three  skid  wires  helically  wound  around  the  outer  surface  of  the  dielectric. 
Preliminary  measurements  indicated  that  the  cable  could  not  be  treated  as  a 
coaxial  cable  with  the  skid  wires  serving  as  the  return  conductor.  The 
calculated  characteristic  impedance  of  a coaxial  cable,  of  the  same  dimensions 
is  about  40  ohms.  The  measured  impedance  varied  between  50  and  75  ohms.  The 
reduced  capacitance. between  the  skid  wires  (as  opposed  to  a solid  coaxial 
sheath)  and  the  center  conductor  results  in  a higher  characteristic  impedance, 

Zg,  since 


ZQ  = kTTC  . 


(1) 
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Another  indication  of  the  non-coaxial  cable  behavior  was  the  unstable 

TDR  results  obtained  on  measurements  of  the  coiled  cable.  The  skid 

wires  did  not  provide  adequate  shielding  between  adjacent  loops  of 

the  cable.  In  order  to  overcome  these  problems  an  improved  outer  conductor 

was  needed  and  a copper  braid  sheath  was  used.  With  the  improved  coaxial -like 

behavior,  the  characteristic  impedance  was  measured  to  be  about  37  ohms  -- 

near  the  predicted  value  of  40  ohms. 

In  a TDR  measurement  a pulse  with  a fast  risetime  is  sent  down  a cable 
and  is  reflected  back  by  variations  in  the  cable  impedance.  The  output  of 
a TDR  measurement  is  the  reflection  coefficient,  p.  A positive  reflection 
referenced  from  the  peak  of  the  incident  pulse  indicates  the  impedance  is 
greater  than  the  TDR  output  impedance.  The  reflection  coefficient  varies 
from  -1  to  1..  The  reflection  coefficient  is  related  to  the  load  impedance 
(R^)  by  the  equation 


rl  = Z0(l  + p)/(l  - p)  . (2) 

This  relation  with  ZQ  equal  to  50  ohms  is  used  in  the  following  discussion 
in  order  to  permit  the  approximate  calculation  of  the  impedance  along  the 
cable  length  from  the  TDR  response. 

Figure  2 shows  the  measured  impedance  (average  over  50  runs)  along 
the  entire  cable  length.  The  response  over  the  first  17  ns  shows  an  impedance 
of  50  ohms.  This  represents  the  impedance  of  the  connector  and  RG-58/U  cable 
from  the  TDR  output  to  the  distribution  cable  input.  The  response  from  about 
17  ns  to  80  ns  shows  the  impedance  of  the  distribution  cable.  Since  the  cable 
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FIGURE  2.  TDR  resDonse  of  full  cable  length. 


-10- 


was  unterminated,  the  impedance  rapidly  increases  at  the  cable  end.  The 
times  shown  are  times  for  the  signal  to  reach  a point  and  return.  Therefore, 
the  time  for  the  TDR  signal  to  propagate  the  length  of  the  cable  is  actually 
about  31.5  ns.  Assuming  a dielectric  constant  of  2.26,  the  length  of  the 
cable  can  be  calculated  to  be  (3  x 10®  m/s)(31.5  x 10  ^s)//2.26  or  6.3  meters, 
close  to  the  actual  length  of  6.2  meters. 

The  circled  region  of  Fig.  2 is  the  point  of  interest  for  the  following 
discussion.  It  represents  a small  segment  of  the  cable  located  about  1.9  meters 
from  the  input  end  of  the  cable.  Figure  3 shows  the  impedance  (average  over 
200  runs)  of  a 0.5  m length  (2.5  ns)  of  the  cable  centered  about  this 
point.  The  impedance  is  seen  to  vary  by  about  0.1  ohm  along  this  distance. 

A mass  of  120  grams  (1.2  newton  force)  had  been  placed  on  the  cable  at 
this  location.  Figure  4 shows  the  measured  impedance  (average  over  200 
runs)  with  and  without  this  force  applied  (lower  and  upper  trace,  respectively). 
Apparently  application  of  the  force  produces  a net  decrease  in  the  separation 
of  the  outer  and  inner  conductors,  a localized  capacitance  increase,  and  a 
consequent  (Eq.  1)  localized  decrease  in  impedance. 

The  difference  between  the  two  traces  is  shown  in  Fig.  5.  The  difference 
has  been  scaled  by  a factor  of  10  to  show  the  sensitivity  available.  The 
circled  region  in  Fig.  5 has  been  further  amplified  (Fig.  6)  in  order  to 
demonstrate  the  ultimate  meaningful  sensitivity  of  this  particular  set  of 
measurements.  The  impedance  variation  over  this  10  cm  length  is  seen  to  be 
about  3 x 10~3  ohm.  Impedance  variations  less  than  2 x 10"^  ohm  (5  ppm)  can 
be  resolved.  With  increased  sampling  frequency  and  improved  experimental 
techniques  an  ultimate  TDR  resolution  of  1 ppm  or  better  in  such  impedance 
measurements  may  be  expected. 
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Figure  3.  TDR  response  of  0.5  m length  cable  segment. 
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FIGURE  4 TOR  response  with  and  without  external  force  (lower  and  upper  trace). 
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Figure  5.  Impedance  difference  from  Fig. 
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FIGURE  6.  TDR  response  of  10  cm  segment  showing  ultimate  resolution. 
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This  series  of  measurements  suggests  that  if  the  dielectric  constant 
changes  by  a part  per  million  over  a few  millimeters  length,  the  TDR  technique 
should  resolve  it. 

The  cable  dielectric  was  purposely  altered  in  order  to  see  if  such 
changes  were  indeed  detectable.  A hole  was  drilled  through  the  cable  sheath 
and  dielectric  to  the  center  conductor.  The  resulting  hole  was  2.3  mm  in 
diameter  and  4.6  mm  deep.  The  TDR  response  of  the  cable  with  the  empty 
hole  was  obtained  and  then  the  hole  was  filled  with  either  isopropyl  alcohol 
(2-propanol)  or  water.  Figures  7 and  8 show  the  impedance  change  resulting 
from  filling  the  hole  with  water  and  alcohol  respectively.  In  Fig.  9 
the  TDR  response  of  a portion  of  the  cable  with  the  empty  hole  (upper-trace) , 
alcohol-filled  hole  (middl e- trace) , and  water-filled  hole  is  shown.  The 
length  of  cable  shown  in'  Figs.  7-9  is  (5  x 10'^s)(3  x 108  m/s)//2. 26  or 

3 

about  1 meter.  The  volume  of  the  dielectric  in  this  cable  is  about  240  cm 
per  meter  of  length.  The  volume  of  the  hole  drilled  through  the  dielectric 

3 

is  about  0.020  cm  . The  resulting  defect  represents  a change  in  the  dielectric 
volume  in  a length  of  1 meter  of  80  ppm.  The  dielectric  properties  of  the 
hole  are  a function  of  the  fluid  it  contains.  Water  has  a static  dielectric 
constant  of  about  80  and  2-propanol,  18.  This  is  consistent  with  the  trend 
shown  in  Fig.  9.  One  difficulty  with  this  experiment  is  that  the  purity  of 
the  liquids  is  questionable  within  this  hole  drilled  in  the  dielectric. 
Particularly  of  concern  is  the  resistivity  of  the  fluid  under  these  conditions. 

This  series  of  measurements  seems  to  imply  that  TDR  techniques  would 
be  effective  in  detecting  at  least  one  type  of  incipient  fault.  In  extruded 
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FIGURE  8.  Impedance  change  caused  by  filling 
cable  dielectric  with  alcohol. 


EXPANDED  LINE  IMPEDANCE 
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FIGURE  9.  TDR  response  of  air-filled  (upper  trace),  alcohol -fi 1 led 
(middle  trace),  and  water-filled  (lower  trace)  2.3  mm 
diameter  hole  in  cable  dielectric. 
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polyethylene  cables  partial  discharge  activity  promotes  the  growth  of  hollow 

channels  in  the  dielectric.  These  channels  branch  out  forming  a tree-like 

structure.  Cables  in  service  for  several  years  have  been  found  to  contain 

3 

many  of  these  "trees".  If  these  "trees"  extend  from  the  inner  to  outer 
conductor,  breakdown  does  occur.  The  inner  walls  of  the  "trees"  are 
conductive^  suggesting  from  the  above  experiment  that  TDR  techniaues  could 
detect  them. 

There  is  a further  difficulty  with  this  series  of  measurements.  It 
appears  in  Fig.  9 that  the  length  of  the  disturbed  area  is  2 nanoseconds 
or  about  0.4  meter  long.  The  actual  hole  diameter  is  only  2.3  mm.  The  water  and 
alcohol  either  must  be  diffusing  a considerable  distance  outside  the 
constraints  of  the  hole,  probably  along  the  stranded  wires  that  constitute 
the  inner  conductor,  or  the  high-frequency  components  of  the  TDR  signal  are 
being  attenuated  so  much  that  spatial  resolution  is  hampered. 

Further  experiments  were  next  initiated  to  determine  if  "trees"  or 
localized  sites  damaged  by  partial  discharges  would  be  discernible  by 
time-domain  reflectometry  techniques.  Measurements  were  performed  on  the 
same  cable  described  above.  A fine  steel  sewing  needle  1.3  mm  in  diameter 
was  carefully  inserted  into  the  dielectric  of  the  cable  to  within  about  one 
millimeter  of  the  center  conductor.  (The  outer  semiconducting  layer  was 
previously  removed  from  the  vicinity  of  the  needle  site).  Negative  pulses 
having  a duration  of  500  ps  were  applied  between  the  needle  and  the  center 
conductor.  At  first,  10  pulses  of  500  volts  amplitude  were  applied  and  the 
TDR  response  obtained.  Then  the  voltage  was  increased  in  500  volt  increments 
and  more  data  were  obtained.  The  maximum  voltage  applied  (the  limit  of  the 
power  supply)  was  2500  volts.  The  field  at  the  tip  of  the  needle  under  these 
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conditions  is  estimated  to  be  of  order  500  kV/mm.  That  is  certainly  high 
enough  to  cause  damage  to  the  dielectric.  However,  after  applying  these 
impulses  even  at  the  maximum  voltage  level,  no  change  was  seen  in  the  TDR 
response. 

While  the  voltage  levels  appeared  to  be  high  enough  to  cause  some 
damage,  the  lack  of  change  in  the  TDR  response  suggested  that  the  dielectric 
v/as  not  being  signif icantly  altered.  In  order  to  test  that  hypothesis,  dc 
voltage  was  applied  to  the  needle.  At  7 kV  a breakdown  was  observed.  The 
TDR  response  showed  no  change.  After  an  additional  10  breakdowns  another 
TDR  response  was  taken.  The  results  are  shown  in  Fig.  10.  The  portion  of 
cable  over  which  the  impedence  is  being  measured  is  (3  x 10^  cm/s)  x 
(2  x 10"  s/2)//2. 26  or  about  20  centimeters  long.  The  reflection  coefficient 
is  negative  because  of  the  cable's  impedance  of  37  ohms  compared  to  the 
output  impedance  of  the  TDR  unit  of  50  ohms.  The  data  shown  in  Fig.  10 
are  the  average  of  512  runs.  Any  change  in  the  TDR  response  is  lost  in  the 
noise. 

The  result  was  initially  surprising.  Previous  work  had  shown  this 
technique  to  be  quite  sensitive.  We  assumed  that  the  partially  conducting 
surfaces  of  the  breakdown  channels  or  "trees"  would  be  readily  discernible  using 
TDR.  To  further  test  the  system,  a hole  1.5  mm  in  diameter  was  drilled  in 
the  cable  dielectric  to  the  center  conductor.  The  hole  was  then  filled 
with  graphite  (pencil  lead).  The  results  are  shown  in  Fig.  11.  The  TDR 
technique  can  readily  detect  the  change  in  impedance.  The  graphite  effectively 
is  forming  a resistive  shunt  between  the  outer  and  inner  conductors.  This  is 
similar  to  the  case  of  a "tree"  extending  across  the  dielectric. 
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FIGURE  1 
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. Change  in  TOR  response  of  cable  resulting  from  11  breakdowns  at 
7 kV  dc.  Breakdown  channel  is  estimated. to  be  1 mm  long  and  is 
undetectable  in  the  noise. 
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FIGURE 
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11.  Change  in  TDR  response  of  cable  resulting  from  filling  a drilled 
hole  in  the  cable  with  graphite.  The  hole  is  1.5  mm  in  diameter 
and  4 mm  deep. 
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The  negative  results  shown  in  Fig.  10  prompted  a study  to  determine 
whether  they  arise  from  a fundamental  limitation  of  the  TDR  or  swept- 
frequency  approach.  In  order  for  a reflection  (or  interaction)  to  occur 
from  a region  of  inhomogeniety  within  a dielectric,  the  minimum  size  of 
that  region  must  be  the  same  order  of  magnitude  as  the  wavelength  of  the 
incident  signal  (e.g.,  an  electromagnetic  signal  with  a minimum  wavelength 
component  of  one  meter  would  have  little  interaction  with  a defect  one 
millimeter  in  length). 

In  the  case  of  our  TDR  measurements,  the  incident  pulse  has  a risetime 
of  less  than  30  picoseconds  implying  a maximum  frequency  component  of  the 
order  of  30  GHz  or  a wavelength  in  the  dielectric  of  about  6 mm.  The  estimated 
length  of  the  breakdown  channel  in  Fig.  10  is  1 mm.  Assuming  the  high 
frequency  components  actually  get  to  the  damage  site,  detection  should  indeed  be 
possible. 

In  order  to  determine  if  the  higher  frequency  components  were  being 
severely  attenuated  making  detection  of  small  sites  impossible,  further 
measurements  were  made.  Figure  12  shows  the  TDR  response  of  a shorted  RG-58/U 
cable  and  the  power  cable  shorted  at  the  damage  site  previously  studied.  Whereas 
the  effect  of  the  short  in  the  RG-58/U  cable  is  evident  over  a time  interval 
of  the  order  of  100  ps  (length  of  1 cm),  the  effect  of  the  short  in  the  power 
cable  extends  over  2 ns  (20  cm  in  length).  The  difficulty  in  resolving  the 
location  of  the  short  in  the  power  cable  demonstrates  that  high  frequency 
attenuation  is  occurring.  Without  these  high  frequency  components,  resolution 
of  small  damage  sites  is  impossible. 

Two  solutions  to  this  problem  have  been  considered:  (1)  better  impedance 

matching  to  increase  the  high  frequency  content  within  the  power  cable;  and 
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TDR  RESPONSE  OF  SHORTED  RG  58  CABLE 
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FIGURE  12.  TDR  response  of  RG-58/U  cable  and  power  terminated  by  zero 
impedance. 
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(2)  after  careful  measurements,  compensation  for  the  attenuation  in  software. 

As  seen  in  Fig.  12  and  discussed  above,  there  is  a serious  attenuation  of  the 
high  frequency  components  of  the  step  input.  A large  fraction  of  this  "loss'1 
may  be  the  result  of  the  impedance  mismatch  between  the  50  ohm  TDR  system 
and  the  37  ohm  high  voltage  cable.  Additionally  at  or  near  the  connection 
to  the  cable  other  mismatches  are  evident  from  the  reflections  present  in 
that  region.  These  results  suggested  the  need  for  an  impedance  matching 
device  to  improve  coupling  of  the  high  frequency  components  into  the  high 
voltage  cable. 

Impedance  transforming  baluns  were  considered,  but  the  idea  discarded 
since  pulse  risetimes  of  25  ns  or  shorter  are  being  considered.  An  investigation 
into  other  impedance  matching  schemes  led  to  the  use  of  tapered  sections 
of  coaxial  transmission  tines  to  achieve  an  impedance  transformation.  A 
design  having  minimal  reflections  over  a wide  bandwidth  employing  a gradual 
exponential  taper  was  examined.  A linear  conical  tapered  line  was  also 
investigated.  Although  the  exponential  taper  is  theoretically  better,  the 
difference  is  so  little  so  as  to  make  the  difficult  exponential  taper  not 

5 

worth  the  extra  machine  time  and  costs.  This  is  especially  true  if  the  taper 
is  long  and  gradual  in  comparison  to  the  wavelength  of  the  lowest  frequency  of 
interest.  The  tapered  section  will  act  as  a high  pass  filter  so  long  as  the 
phase  constant  does  not  contain  an  imaginary  term.  At  lower  frequencies  the 
phase  constant  is  complex  and  the  signal  is  attenuated.  Practically, 
impedance  changes  of  2 to  1 can  be  made  if  the  tapered  section  is  one  wavelength 
long  or  longer.  A good  impedance  match  will  result  for  higher  frequencies. 
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Figure  13  shows  details  of  the  tapered  line.  The  outer  conductor 
of  the  line  is  made  from  a thick-walled  section  of  brass  tubing 
which  clamps  over  the  braided  shield  of  the  high  voltage  cable. 

In  the  right  end  of  the  "tube",  a machined  adaptor  is  inserted.  This 
couples  to  a General  Radio^  type  874-BP58A  flanged  panel  connector.  This 
type  of  connector  is  known  to  have  a constant  50-ohm  impedance  throughout 
its  length  and  to  provide  for  a good  electrical  connection  with  low  series 
and  shunt  losses.  The  stepped  portion  of  the  center  conductor  shov/n  in 
Fig.  13  fits  inside  the  adaptor.  The  dimensions  have  been  selected  to 
maintain  a 50-ohm  characteristic  impedance  throughout  its  entire  length. 

The  susceptances  at  each  "step"  should  be  low  enough  even  at  the  highest 
frequencies  to  provide  a constant  50-ohm  match  in  this  region  of  the  line. 

The  center  conductor  of  the  line  is  tapered.  The  tapered  length  is 
30.18  cm  which  is  one  wavelength  at  a frequency  of  993  MHz.  It  would  be 
desirable  to  have  a longer  line  so  the  lower  cutoff  frequency  would  be  less, 
but  this  becomes  mechanically  difficult  and  cumbersome. 

The  left  end  of  the  tapered  center  conductor  fits  directly  over  the 
center  conductor  of  the  high  voltage  cable.  Because  of  possible  end  effects 
of  the  tapered  section  and  the  abrupt  discontinuity  of  the  Dolyethylene 
dielectric  at  the  end  of  the  cable,  impedance  mismatch  reflections  can  be 
expected.  The  influence  of  these  remaining  reflections  has  been  experimentally 
determined. 


50-OHM  COAXIAL 
CONNECTOR  - 
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OVER  OUTSIDE  OF  CABLE 
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The  performance  of  the  transition  section  is  shown  in  Figs.  14  to  17. 

Figure  14  shows  the  TDR  response  at  the  end  of  the  50-ohm  coaxial  cable  and  at 
the  beginning  of  the  37-ohm  distribution  cable  without  the  transition  section. 
The  output  of  the  TDR  analyzer  was  adapted  to  BNC-type  coaxial  connectors. 

A one  meter  length  of  RG-58/U  cable  was  connected  to  a GR^  type  874  BNC-to-GR 
adaptor.  The  six  outer  strands  of  the  distribution  cable  center  conductor 
were  removed  leaving  only  one  strand  at  the  center.  The  female  portion  of 
the  BNC-to-GR  adaptor  was  placed  over  the  single  remaining  center  strand  of 
the  cable's  inner  conductor.  The  braided  sheath  of  the  cable  was  then  dressed 
over  the  adaptor  and  the  whole  assemblage  taped  into  place.  Using  this 
arrangement,  two  major  impedance  mismatches  occurred.  The  first  was  due 
to  the  short  length  of  small  diameter  single-strand  center  conductor.  This 
length  was  about  9 mm  and  caused  the  impedance  to  increase  to  about  90  ohms, 
the  first  major  peak  in  Fig.  14.  The  second  mismatch  is  due  to  the  lower 
characteristic  impedance  of  the  distribution  cable.  (The  small  variations 
of  impedance  seen  in  the  37-ohm  portion  of  the  cable  are  caused  principally 
by  the  unevenness  of  the  braided  sheath). 

Figure  15  shows  a cross-sectional  view  of  the  50-ohm  end  of  the 
transition.  (The  GR874BP  connector  is  not  shown  on  the  left  end).  The 
steps  of  the  center  conductor  were  designed  to  coincide  with  the  steps  of 
the  outer  conductor.  However,  an  error  on  the  design  drawing  of  0.2  cm  shifted 
the  two  larger  diameters  on  the  right  towards  the  left  end  resulting  in  a 
lack  of  coincidence  between  the  inner  and  outer  conductors.  This  was  discovered 
after  the  piece  was  fabricated  and  assembled  and  the  TDR  response  obtained. 
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gure  14.  TDR  response  of  RG-58/U  cable  connection  onto  37-ohm  distribution 
cable. 
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Figure  16.  50  ohm  end  of  the  transition  section  per  second  design. 
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The  TDR  response  is  keyed  to  the  stepped  center  conductor  as  shown 
in  Fig.  15.  (Vertical  sensitivity  is  about  5 ohms  per  division).  Several 
problems  were  evident  upon  examination  of  this  response.  First  the  impedance 
variations  were  greater  than  expected  and  second,  the  impedance  of  the 
fourth  step  was  lower  than  design  specifications.  From  the  latter,  a 
second  drawing  error  was  discovered:  the  inside  diameter  of  the  outer 

conductor  was  too  small  by  0.09  cm. 

It  was  obvious  that  even  though  the  steps  in  diameter  were  each  only 
0.12  cm,  a 3-to  5-ohm  impedance  change  occurred  at  each  step.  Because 
the  impedance  decreased  in  the  vicinity  of  each  step,  this  indicated  that 
the  shunt  susceptance  at  the  steps  was  too  great.  In  order  to  compensate 
for  and  to  reduce  this  susceptance,  each  of  the  steps  of  the  center  conductor 
piece  was  shifted  to  the  "right  by  an  amount  equal  to  D/8  where  D is  the 
inside  diameter  of  the  outer  conductor.  Also  the  inside  diameter  of  the  last 
step  of  the  outer  conductor  was  increased  by  0.09  cm  to  increase  the  impedance 
in  that  region  of  the  transition  section.  The  TDR  response  in  Fig.  16  shows 
a significant  improvement  over  that  of  Fig.  15.  The  impedance  does  not  vary 
by  more  than  two  ohms. 

Figure  17  shows  the  TDR  response  for  the  entire  length  of  the  tapered 
transition  section.  At  the  left  is  the  50-ohm  TDR  system  impedance;  in  the 
center  is  the  tapered  line;  and  on  the  right  is  the  37-ohm  high  voltage 
distribution  cable.  As  can  be  seen  at  the  right  end  of  the  tapered  section, 
impedance  variatiations  are  apparent.  These  are  caused  by  the  end  effects  of 
the  tapered  line  attaching  to  the  distribution  cable  center  conductor  and 


Figure  17. 


TDR  response  of  tapered  line  transition  section. 
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imperfections  in  the  shield  geometry  at  and  near  the  beginning  of  the 
distribution  cable. 

Comparing  Figs.  14  and  17  we  see  that  the  transition  section  has 
fulfilled  its  intended  purpose  --  a smooth  impedance  transition  from  the 
50-ohm  output  impedance  of  the  TDR  analyzer  to  the  37-ohm  characteristic 
impedance  of  the  cable.  The  next  step  was  to  determine  if  this  transition 
section  permitted  significantly  more  high  frequency  signal  to  enter  the 
cable  and,  if  so,  would  this  enable  resolution  of  incipient  faults  or 
would  attenuation  within  the  cable  still  be  a problem. 

With  the  transition  section  in  place  the  outer  and  inner  conductors 
were  shorted  together  at  a distance  of  61  cm  from  the  cable  end  (and  the 
TDR  analyzer).  The  TDR  display  was  averaged  512  times  and  sent  to  the 
intelligent  terminal  for  'display.  The  results  are  shown  in  Fig.  18.  The 
ordinate  has  not  been  normalized  as  on  the  previous  figures  (i.e.,  a 
constant  has  not  been  added  to  the  reflection  coefficient  in  order  to 
force  it  to  be  zero  for  the  power  cable's  characteristic  impedance  of  37  ohms) 
The  important  parameter,  reflection  coefficient  variation,  is  not  affected. 

The  upper  graph  is  from  Fig.  12  and  is  taken  without  benefit  of  the  transition 
section.  The  lower  graph  demonstrates  the  effect  of  the  transition  section. 
For  the  sake  of  simplicity  we  will  define  the  decay  time  as  the  time  it  takes 
the  signal  to  decrease  from  7/8  to  1/8  of  the  difference  of  the  average 
maximum  and  minimum  values.  The  decay  time,  thus  defined,  is  1750  ps  without 
and  800  ps  with  the  transition  section.  This  means  the  additional  high 
frequency  components  entering  the  cable  with  the  transition  section  have 
permitted  better  spatial  resolution  of  the  short.  The  800  ps  decay  time 
corresponds  to  a length  of  8 cm  which  is  much  longer  than  the  shorted  hole's 
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' figure  18.  jqr  response  of  power  cable  terminated  by  zero  impedance  with  and 
! without  transition  section. 
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diameter  (0.13  cm).  Sufficient  high  frequency  content,  therefore,  is 
still  not  available  to  resolve  this  region  adequately. 

The  second  factor  affecting  high  frequency  propagation  is  the  attenuation 
in  the  cable  itself.  Figure  19  shows  results,  taken  as  before,  with  the  short 
located  at  two  different  distances  from  the  TDR  analyzer.  The  figure  indicates 
that  the  greater  the  distance  of  propagation  in  the  cable  the  greater  the 
decay  time  of  the  TDR  response.  A short  located  at  4 cm  from  the  cable  end 
has  a decay  time  of  340  ps  while  a short  located  at  163  cm  has  a decay  time 
of  1050  ps.  From  above  we  found  that  a short  61  cm  from  the  cable  end  has 
a decay  time  of  800  ps.  At  4 cm  the  decay  time  corresponds  to  a spatial 
resolution  of  3.4  cm  and  at  163  cm  to  10.5  cm.  These  results  are  troubling 
in  that  they  imply  that  the  intrinsic  high  frequency  attenuation  in  this 
type  of  cable  makes  the  detection  at  large  distances  of  small-size  damage 
sites  very  difficult,  if  not  impossible. 

Figure  20  shows  this  quite  clearly.  In  the  previous  measurements,  we 
had  not  been  able  to  detect  a hole  drilled  through  the  dielectric  1.3  mm  in 
diameter.  We  again  attempted  this  with  holes  drilled  4 cm  from  the  cable 
end  and  61  cm  from  the  cable  end.  At  61  cm  the  response  with  the  hole 
(light  trace)  and  without  (heavy  trace)  are  nearly  the  same,  any  differences 
not  being  repeatable.  At  4 cm,  there  is  a definite  difference  between  the 
two  traces.  New  minima  and  maxima  have  been  introduced  with  the  drilling  of 
the  hole.  From  these  data,  it  is  inferred  that,  at  distances  greater  than 
a few  centimeters,  the  high  frequency  attenuation  has  caused  this  information 


to  be  lost. 
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SHORTED  POWER  CABLE,  LENGTH  4cm 


TIMECSECONDS3 

SHORTED  POWER  CABLE,  LENGTH  163cm 


TIHECSECONDSD 

. 

j Figure  19.  TDR  response  of  power  caDle  terminated  by  zt>'n  impedance;  the 

location  of  short  is  4 cm  and  163  cm  from  cable  end;  measurements 
made  using  transition  section. 


REFLECTION  COEFFICIENT  REFLECTION  COEFFICIENT 


TDR  RESPONSE  OF  CABLE  AT  Aem 


TIME  CSECCNDSD 

TDR  RESPONSE  OF  CABLE  AT  61cm 


0.0080 


0.0000 


-0.0080 


0 . 0E+000 


S . 0E-0 1 0 
TIME  CSEC0NDS3 


1 . 0E-003 


TOP  response  of  power  cable  (using  transition 
(lighter)  and  without  (darker)  a hole  drilled 
dielectric.  Distance  of  1.3  mm  diameter  hole 
is  4 cm  (upper)  and  61  cm  (lower). 


section)  with 
through 

from  cable  end 


Figure  20. 
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The  results  indicate  that  incipient  fault  detection  schemes  based  on 
high  frequency  sounding  may  not  be  practical  for  plastic-insulated  cables 
although  it  is  possible  that  such  techniques  may  be  of  value  for  compressed-gas 
(spacer  failure),  superconducting,  and  other  cable  types.  It  is  also  possible 
that  our  results  obtained  on  a 15  kV  extruded  polyethylene  distribution  cable 
may  not  fully  apply  to  higher  voltage  transmission  cables  of  similar  design. 

High  frequency  sounding  could  still  be  useful  for  fault  location  and,  if  a 
sufficiently  large  damage  site  is  available,  incipient  fault  location. 

1 1. 3 Software  Development 

In  addition  to  the  software  presented  in  Appendix  A which  was  developed 
for  data  transfer  between  the  digital  processing  oscilloscope  and  the 
intelligent  terminal,  other  software  needed  to  be  developed.  Experiments 
in  the  coming  year  will  be  initiated  with  the  objective  of  determining  the 
signal  velocity  and  attenuation  in  cables  at  rf  frequencies.  Initial 
measurements  will  be  in  the  time-domain.  Step-like  voltage  waveforms 
will  be  applied  to  the  cable  and  the  output  waveform  digitized  and 
transmitted  to  the  smart  terminal.  A Fourier  transform  will  then 
yield  the  transfer  function  in  the  frequency  domain.  Problems  arise, 
however,  because  instead  of  a Fourier  transform  of  continuous  data  we 
have  a Fourier  transform  of  discretely  sampled  data  or  a discrete  Fourier 
transform  (DFT).  A DFT  of  an  impulse-like  waveform  (a  waveform  which  ends 
at  the  same  amplitude  it  begins)  is  straightforward  and  a fast  algorithm 
exists  for  its  solution  (FFT).^  Step-like  waveforms,  however,  will  produce 
a serious  truncation  error'7  if  they  are  treated  by  conventional  DFT  (and  FFT) 
techniques. 

We  have  found  two  different  (but  related)  methods  to  overcome  this  problem. 

8 

The  first  is  the  sin  x/x  expansion  presented  by  Samulon. 
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The  measured  curve  is  sampled  at  time  intervals,  x,  such  that: 


x<l/2fc  . 


(3) 


where  f is  the  highest  frequency  component  present.  This  restraint  is 

7 9 

necessary  to  avoid  "aliasing"  problems  ’ . Using  a sin  x/x  expansion  the 
given  time  function,  F(t),  can  be  expressed 


00  A sin  f (t-nx)2ir 
F(t)  = A '--f  it-nx)2, 


(4) 


n=o  c 


where  An  is  the  amplitude  of  the  nth  sampled  point.  The  Fourier  transform 
of  the  nth  term  of  Eq.  (4)  is  given  by 


.<->  -I  ^ 


A sin[2rrf  (t-nx)]  e'^^dt 


2 Trf  c ( t-nx ) 


(5) 


or 


~ 2tt?c J 


A^  r sin[2Trf  (t-nx)]  {cos  wt  - j sin  cot}  dt 


t-nx 


(6) 


Subsituting  t'  for  t-nx  and  using 


ex 

/ 


sin  mx  cos  nx  dx 


= it  for  m>n>o 


(7) 
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We  obtain 


4>n(w)  - 2f 

c 


(8) 


or 


c 


(9) 


Software  based  or.  Eq.  (9)  has  been  written  and  is  incorporated  in 
Appendix  B.  However,  this  is  not  the  most  useful  form  of  the  sin  x/x  expansion. 
If  An  does  not  equal  zero  for  some  n greater  than  an  arbitrary  m,  Eq.  (9)  does 
not  converge  and  cannot  be  used.  If  An  does  not  converge  and  the  summation  is 
arbitrarily  terminated  at  some  point,  ripples  will  occur  obscuring  the  results. 

A modified  form  of  Eq.  (9)  can  be  used  which  does  not  have  such  limitations. 
Let  Bn  = An  - An_.j  and  replace  the  sequence  {AQ,  A-j , A^  ...}  by  {AQ,  A-j  e~T//<5 , 

. . }.  The  value  of  5 can  be  made  arbitrarily  large  so  as  not  to  affect 
our  transient  unless  i is  very  large  or  in  the  frequency  domain  oo  is  very  small. 
Using  the  relation 


n 


(10) 


Eq.  (9)  becomes 
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$(u)) 


E 

n=o 


n 

( E 


m=o 


B )e“JnwT  “ nT/5 
m' 


(ID 


If  we  let  ne  = jnoox  + nx/6,  we  obtain 


2f  $(uj) 
c 


00  n 

E ( E Bm)e' " 
m 

n=o  m=o 


(12) 


or  equivalently 


oo 

00 

2f  Ha)  = ( E e‘ne)  ( E B e'mE) 
n=o  m=o 


(13) 


~ - n c 

The  series  Ee  is  absolutely  convergent  if 


^ 1 0- jnarr  - naj/5  . 


(14) 


converges.  But 


e-jnu)T  - nx/6 , = g-nx/6 


(15) 


and  (14)  is  seen  to  be  a geometric  progression,  convergent  for  e~T/  <1 . 
Using  the  formula  for  the  sum  of  an  infinite  geometric  progression: 


CO 


E e 
n=o 


= 1/0  - e"e) 


(16) 
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o r 

l e"ne  = 1/(1  - e‘T/5e'J“T)  . (17) 

n=o 

The  factor  e~T//(5  can  be  made  arbitrarily  close  to  unity  yielding 

E e"ne  = 1/(1  - e"ja3T)  . (18) 

n=o 

From  Eq.  (13),  we  obtain 

. oo 

Urn)  = [l/(2f  (1  - e'J“T)]  2 Bn  e'0™1  . (19) 

n=o 

Appendix  C contains  the  software  utilizing  Eg.  (19).  This  program  is  nood  for 
waveforms  of  arbitrary  shape  and  can  always  be  used  assuming  t £ l/2f  and  that, 
for  n greater  than  some  m,  Bn  approaches  some  constant  (not  necessarily  zero). 

The  only  difficulty  with  Eg.  (19)  is  that  the  processing  can  take  several  seconds 
per  frequency  point  making  it  impractical  for  some  experiments. 

The  second  method  used  a modified  FFT.  As  mentioned  above  a conventional 
FFT  yields  incorrect  results  due  to  truncation  error.  A modified  FFT  has 
recently  been  presented^  which  overcomes  this  problem. 

From  Eq.  (19)  and  the  definition  of  Bn  we  obtain 

tU)  * [l/(2fcO  - e'jVr))]  z (A  - A )e'jn“  . (20) 

n=o 

We  assume  the  waveform  becomes  a constant  at  T = Nr  (i.e.,  A..  = A.  ^ for  i>N). 

It  is  further  assumed  that  AQ  is  equal  to  zero  (i.e.,  the  dc  offset  of  the 
waveform  is  zero.  A . is  defined  to  be  zero.  Substitutinq  v 

-I r 


for  v 
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(v  = co/2n)  such  that 


v = kv  = ic(l/T)  = k/Nt,  k = 1,2,  N - 1 
< o 


Eq.  (20)  becomes 


<±> 

K 


[t/(1  - e-jK2¥/N)]V  (A  - a )e-jm<*T/N 

n=o  n n 


(21) 


(22) 


In  Eq.  (22),  t has  been  set  equal  to  l/2f  satisfying  the  sampling  theorem. 
In  order  to  be  consistent  with  FFT  conventional  formalism  the  factor  x will 
be  removed  from  Eq.  (22).  Separating  the  sum  into  two  parts,  we  obtain 


= 1/(1  - e~Jl<2Tr/N)  {V  A„e-jnK2,T/N  - V An_ie-J'nK2,r/N}  . (23) 


n=o 


n=o 


Defining  S , 


S 


K 


Z A e-jnK2,r/H  , 


(24) 


to  be  the  FFT  acting  on  the  data  set  { An } , Eq.  (23)  becomes 

$ - i/(i  . e-jx2VN}  {s  . N’1  A e-jn<2ir/N}  , (25) 

k k n-1 


or 


i,  = (1  - e-j^^1)-1  (S 

< ' K 


e~  j k2tt/N 


N-2 

Z 

9=0 


V 


- j J?,k2tt/N 

} 


(26) 
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This  is  equivalent  to 


= O - e 


and,  since  * 1, 


4-  = S + Am  ,/(l  - e 
k k N-r 


-jk2tt/N 


) k = 1 , . . . , M - 1 . 


(28) 


This  is  a very  important  result  in  that  it  permits  considerable  saving 
in  computational  time.  Using  an  FFT  read-only-memory  pack  our  intelligent 
terminal  can  compute  for  256  frequency  points,  in  under  a minute  compared 
to  over  an  hour  if  Eq.  (19)  were  used.  Appendix  D contains  the  software 
utilizing  Eq.  (28).  Appendix  E contains  the  results  of  using  the  sin  x/x 
approach  Eq.  (19)  and  the  modified  FFT  approach  Eq.  (28)  on  the  same  data.  The 
results  are  identical. 

Both  approaches  are  sensitive  to  dc  offset.  Assuming  the  zero 
frequency  result  is  not  of  interest  the  waveform  should  be  shifted  so  as 
to  make  the  amplitude  zero  at  time  equal  to  zero.  Equivalently,  using  Eq.  (19), 
Bq  can  be  set  to  zero  in  software  or  using  the  modified  FFT  approach  replacing 
Eq.  (28)  by 


(29) 
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If  the  software  is  modified  in  such  a way  both  methods  are  insensitive,  as 
they  should  be,  to  dc  offset. 

The  software  development  will  be  invaluable  in  future  measurements. 

The  sin  x/x  approach  has  the  advantage  of  allowing  the  user  to  choose  the 
points  to  be  calculated.  The  modified  FFT  approach  is  considerably  faster. 

III.  Conclusion 

A measurement  program  has  been  initiated  at  NBS  that  will  contribute 
to  the  development,  testing,  and  utilization  of  incipient  fault  detection/ 
location  instrumentation.  The  measurement  system,  centered  on  an  intelligent 
terminal,  permits  the  digitizing  of  transient  and  continuous  waveforms  and 
subsequent  processing. 

The  system  was  initially  configured  in  such  a way  as  to  permit  time-domain 
reflectometry  studies  on  -a  15  kV  distribution  cable.  The  results  were 
significant  in  that  they  showed  that  TDR  techniques  would  have  very  limited 
use  for  detection  of  incipient  faults.  The  difficulty  is  that  the  high 
frequency  attenuation  in  the  extruded  polyethylene  insulated  distribution  cable 
prevents  the  spatial  resolution  of  the  anticipated  incipient  faults. 

Measurements  in  the  future  will  determine  the  rf  properties  of  high 
voltage  transmission  cables.  These  results  will  be  important  in  the  design 
of  incipient-fault  intrumentation . These  measurements  will  be  made  in  the 
time  domain  requiring  mathematical  transformation  to  the  desired  frequency 
domain.  Extensive  software  has  been  written  and  tested  using  two  different 


methods . 
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APPENDIX  A 


This  software  permits  the  collection  of  data  from  the  digital  processing 
oscilloscope  DPO.  The  data  can  first  be  averaged  up  to  4096  times  by  the 
DPO.  The  data  are  plotted  on  the  terminal  screen.  Portions  may  be  expanded. 
The  data  can  be  filtered  (adjacent  points  averaged).  The  data  can  be 
output  to  a plotter  for  quality  graphics. 
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APPENDIX  B 


This  software,  using  the  sin  x/x  approach  described  in  the  text,  Eq. 
can  be  used  for  a discrete  Fourier  transform  of  time-doinain  results. 


(9), 

It  is 


essential  that  the  last  digitized  point  of  the  waveform  equal  zero  or 
truncation  errors  will  result. 
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APPENDIX  C 


This  software  using  a modified  sin  x/x  approach,  Eq.  (19),  results  in  a 
discrete  Fourier  transform  of  time-domain  data.  The  final  point  of 
the  digitized  waveform  need  not  equal  zero  but  the  waveform  is  required  to 
approach  some  constant  value. 
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APPENDIX  D 

This  software  performs  a fast  Fourier  transform  of  time 
technique  is  a modified' Eq.  (28)  Cool ey-Tukey  technique 
transform  waveforms  that  do  not  have  the  same  beginning 
The  waveform  must  begin  at  an  amplitude  of  zero. 
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APPENDIX  E 

Five  different  sets  of  data  are  included.  The  first  two  show  the  results 
of  using  the  software  contained  in  Appendices  B and  C on  the  same  data. 
The  results  are  identical  because  the  waveform  ends  at  zero  amplitude. 

The  other  three  sets  of  data  illustrate  the  performance  of  the  different 
software  approaches,  Appendices  B,  C,  and  D,  to  the  rising  step  (starting 
at  zero,  ending  above  zero).  Only  the  software  contained  in  Appendices  C 
and  D handle  the  data  correctly. 
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The  following  three  pages  demonstrate  the  results  of  applying  the  software 
contained  in  Appendix  B to  a step  with  endpoint  equal  to  zero. 

The  first  page  shows  the  waveform  as  acquired  from  the  DPO  and  a reduced 
form  of  the  waveform  which  includes  only  every  fourth  data  point.  The 
sin  x/x  software  yields  the  frequency  spectrum  of  the  reduced  waveform 
on  the  second  page  and  a plot  of  this  spectrum  is  on  the  third  page. 

This  software  correctly  transforms  waveforms  of  this  type. 
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The  following  three  pages  demonstrate  the  results  of  applying  the  software 
contained  in  Appendix  C to  a step  with  endpoint  equal  to  zero. 

The  first  page  shows  the  waveform  as  acquired  from  the  DPO  and  a waveform 
showing  the  difference  between  every  fourth  data  point.  The  modified 
sin  x/x  software  yields  the  frequency  spectrum  of  the  reduced  waveform 
on  the  second  page  and  a plot  of  this  spectrum  is  on  the  third  page. 

This  software  correctly  transforms  waveforms  of  this  type. 
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The  following  ten  pages  show  the  results  of  applying  the  software  contained 
in  Appendix  B to  a step-like  waveform  with  endpoint  not  equal  to  zero. 

The  waveform,  as  acquired  from  the  DPO,  is  on  the  first  page.  The 
transformed  spectra  is  tabulated  on  the  next  eight  pages  and  plotted 
on  the  last  page. 

From  theoretical  considerations  in  the  text  and  results  shown  next  in 

this  appendix,  this  software  is  clearly  not  applicable  to  waveforms  of  this 

type. 
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The  following  ten  pages  show  the  results  of  applying  the  software  contained 
in  Appendix  C to  a step-like  waveform  with  endpoint  not  equal  to  zero. 

The  first  page  shows  the  waveform  as  acquired  from  the  DPO  and  the  difference 
waveform  obtained  by  subtracting  successive  points  of  the  original  waveform. 
The  frequency  spectrum  of  this  difference  waveform  is  tabulated  on  the  next 
eight  pages  and  plotted  on  the  last. 


This  software  correctly  transforms  waveforms  of  this  type. 
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The  following  ten  pages  show  the  results  of  applying  the  software  contained 
in  Appendix  D to  a step-like  waveform  that  begins  at  zero  and  ends  at  some 
nonzero  value. 

The  waveform,  as  acquired  by  the  DPO,  is  shown  on  the  first  page.  The 
transformed  spectrum  is  tabulated  on  the  next  eight  pages  and  plotted  on 
the  last. 


This  software  correctly  transforms  waveforms  of  this  type. 
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